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Abstract 

We model the dark sector of the cosmic substratum by a viscous fluid with an equation of state 
p = — C©) where is the fluid-expansion scalar and C, is the coefficient of bulk viscosity for which we 
assume a dependence C, en on the energy density p. The homogeneous and isotropic background 
dynamics coincides with that of a generalized Chaplygin gas with equation of state p = —A/p"-. 
The perturbation dynamics of the viscous model, however, is intrinsically non-adiabatic and qual- 
itatively different from the Chaplygin-gas case. In particular, it avoids short-scale instabilities 
and/or oscillations which apparently have ruled out unified models of the Chaplygin-gas type. We 
calculate the matter power spectrum and demonstrate that the non-adiabatic model is compatible 
with the data from the 2dFGRS and the SDSS surveys. A x^-analysis shows, that for certain 
parameter combinations the viscous-dark-fluid (VDF) model is well competitive with the ACDM 
model. These results indicate that non-adiabatic unifled models can be seen as potential contenders 
for a General- Relativity-based description of the cosmic substratum. 
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I. INTRODUCTION 



Most current cosmological models rely on the assumption that the dynamics of the Uni- 
verse is described by Einstein's General Relativity (GR) and a material content that is 
dominated by two so far unknown components, pressureless dark matter (DM) and dark 



energy (DE), a su 
tual situation see 



3stance equipped with a large negative pressure. For reviews of the ac- 
2I, l3| and references therein. The preferred model is the ACDM model 
which also plays the role of a reference model for alternative approaches to the DE problem. 
While the ACDM model can describe most of the observations, there still remain puzzles 
{4]. Other attempts to describe the apparently observed accelerated expansion of the Uni- 
verse (see, however, jsl) are a potential back- reaction mechanism from non-linear structure 
formation 61 or models of the Lemaitre-Tolman-Bondi type which are isotropic but not ho- 
mogeneous 7|. Then there exists a line of investigation that modifies GR with the aim to 
obtain an accelerated expansion of the Universe as a result of the (modified) geometrical sec- 
tor instead of matter with a negative pressure [8|, |9|. Here we focus on a class of approaches 
within GR that do not separate DM and DE from the start but regard the dark sector as a 
one-component substratum which exhibits properties of both DE and DM on a joint foot- 
ing. The best known models of this kind are Chaplygin-gas t ype cosmologies. Starting with 



lOj], there has been a considerable activity in this field 11 



12 



Q,Q,[ii 



16|, 



17 



18 



Among the host of models proposed over the years for the dark sector, the unified models are 
minimal in the sense that they assume just one component that describes both dark matter 
and dark energy. Both components manifest themselves observationally only through their 
gravitational action. Therefore, a unified description is certainly attractive, at least as long 
there is no direct, i.e., other than gravitational, detection of either or even both of the dark 
components. While the Chaplygin-gas models could well describe the SNIa results [20|, i.e., 
the cosmic background dynamics, they seem to have fallen out of favor because of appar- 
ent problems to reproduce the matter power spectrum. The difficulties of the generalized 
Chaplygin-gas cosmologies are related to the values of the sound speed of these models. 
Depending on the a-parameter, the sound speed is either finite, i.e., it becomes of the order 
of the speed of light, or its square is negative. In the first case, the small-scale perturbation 
behavior is oscillatory, in the second case there appear instabilities. In neither of these cases, 
the observed matter power spectrum is reproduced. This circumstance has led the authors 



of 



2l| to the conclusion that Chaplygin-gas models of the cosmic medium are ruled out as 



competitive candidates. Similar results were obtained from the analysis of the anisotropy 
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spectrum of the cosmic microwave background [22|, l23|, except possibly for low values of 
the Hubble parameter 2J]. However, these conclusions rely on the assumption of an adi- 
abatic cosmic medium. It has been argued that there might exist entropy perturbations, 
so far not taken into account, which may change the result of the adiabatic perturbation 
analysis 25|, |26|. A problem here is the origin of non-adiabatic perturbations which should 



reflect the internal structure of the cosmic medium. The latter is unknown but it may well 
be more complicated then suggested by the usually applied simple (adiabatic) equations of 
state. Non-adiabatic perturbations will modify the adiabatic sound speed. The speed of 
sound has generally attracted interest as a tool to discriminate between different dark en- 



ergy models 
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3l| . "Silent" Chaplygin gases were postulated by introducing ad 



IOC a non-adiabatic counter-term to exactly compensate the adiabatic pressure contribution 



26 



32|. 



Another option for a unified description of the dark sector are viscous models. From early 
universe cosmology it is known that a bulk viscosity of the cosmic medium can induce an 
inflationary phase [33|. A bulk viscous pressure in the early universe can be the result of 



cosmological particle production 



34 



35|. Under the conditions of spatial homogeneity and 



isotropy, a scalar bulk viscous pressure is the only admissible non-equilibrium phenomenon. 
The cosmological relevance of bulk viscous media has subsequently been investigated in 



some detail for an inflationary phase in the early universe (see 36 



37 



38 



391] and references 



therein). With accumulating evidence for our present Universe to be in a stage of accelerated 
expansion, an effective bulk viscous pressure was discussed as one of the potential sources 



for this phenomenon as well. It was argued in [25|, |40|, that such a pressure can play 
the role of an agent that drives the present acceleration of the Universe. The option of 
a viscosity-dominated late epoch of the Universe with accelerated expansion was already 
mentioned in 4l|] , long before the direct observational evidence through the SN la data. For 
a homogeneous and isotropic universe, the ACDM model and the (generalized) Chaplygin- 
gas models can be reproduced as special cases of this imperfect fluid description [25|. The 
possibility of using cosmological observations to probe and constrain imperfect dark-energy 
fluids was investigated in 42| and 43 1. 

It is obvious that the bulk viscosity contributes with a negative term to the total pressure 
and hence a dissipative fluid seems to be a potential dark energy candidate (For a recent 
preprint see 4J]). However, it is expedient to repeat a cautionary remark. In traditional 
non-equilibrium thermodynamics the viscous pressure represents a (small) correction to the 
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(positive) equilibrium pressure. This is true both for the Eckart 45| and for the Israel- 



Stewart theories 



46 



471]. Here we shall admit the viscous pressure to be the dominating 



part of the pressure. This is clearly beyond the established range of validity of conventional 
non-equilibrium thermodynamics. Non-standard interactions are required to support such 
type of approach 25|, |40|. Of course, this reflects the circumstance that dark energy is 
anything but a "standard" fluid. There are suggestions that viscosity might have its origin 
in string landscape [48|. To successfully describe the transition to a phase of accelerated 
expansion, preceded by a phase of decelerated expansion in which structures can form, it 
is necessary that the viscous pressure is negligible at high redshifts but becomes dominant 
later on. 

Extending previous work 



49 



50j), this paper provides a detailed study of perturba- 
tions in a viscous fluid model of the dark sector which in the homogeneous and isotropic 
background coincides with the dynamics of a generalized Chaplygin gas. The differences in 
the perturbation dynamics of both approaches are traced back to an inherent non-adiabatic 
behavior of the viscous model. In a sense, the viscous dark-fluid model can be seen as a non- 
adiabatic Chaplygin gas. It does not suffer from the shortcomings of the so far considered 
adiabatic Chaplygin-gas models, i.e., it predicts neither (unobserved) oscillations nor insta- 
bilities. Moreover, the non-adiabatic behavior is part of the model and no ad hoc introduced 
counter terms to the adiabatic sound speed are required. The resulting power spectrum fits 
both the 2dFGRS as well as the SDSS data and, for certain parameter combinations, the 
X^- value of the viscous model is better than the corresponding ACDM-value. 

As in Q, [5^, we shall describe the bulk viscous pressure by Eckart 's expression 45] 
p = —S,u\ (Latin indices run from to 3), where the (non-negative) quantity ^ is the 
(generally not constant) bulk-viscosity coefficient and is the fluid-expansion scalar which 
in the homogeneous and isotropic background reduces to SiY, where if = ^ is the Hubble 
parameter and a is the scale factor of the Robertson- Walker metric. By this assumption 
we ignore all the problems inherent in Eckart 's approach which have been discussed and 



resolved within the Israel-Stewart theory 



46|, |47| (see also 



36 



37 



38 



39] and references 



therein). We expect that for the applications we have in mind here, the differences are of 
minor importance. 

The paper is organized as follows. Section [ITl introduces the basic dynamics of bulk 
viscous cosmology and describes the homogeneous and isotropic background solutions in 
analogy to their well-known (generalized) Chaplygin-gas counterparts. Section IIIII contains 
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the general non-adiabatic perturbation dynamics. In order to point out the differences 
between the viscous and the Chaplygin-gas models, the dynamics is developed in parallel 
for both approaches. In Section IIVI it is shown that the basic second-order equations only 
coincide deep in the matter-dominated era. Subsequently, the numerical solutions and the 
corresponding matter power spectra are presented and compared with the 2dFGRS and the 
SDSS data. A summary of our study is given in Section |Vl 



II. BASIC DYNAMICS 



We assume the cosmic medium to be described by an energy-momentum tensor 



T 



ik 



(3) 



with the equation of state 

P = -ce (2) 

of a bulk- viscous fluid, where G = u^.^ is the expansion scalar and ( is the coefficient of bulk 
viscosity. In the homogeneous and isotropic background one has G = 3H, where H is the 
Hubble rate. If, moreover, the background is spatially flat, the Friedmann equation 

3H^ = SnGp 

implies G oc p^/^. Let us further assume that C, (x p" . This corresponds to a background 
equation of state 

with a constant A > 0. Comparing this with the equation of state of a generalized Chaplygin 
gas (subscript c), 

Pc = -^, (5) 

the correspondence a = — (z/ + |) is obvious. This will allow us to apply known results 
from the Chaplygin-gas dynamics to the background dynamics of the viscous model. The 



(4) 



similarity between generalized Chaplygin gases and bulk- viscous fluids is well known |50l. |51|. 
We recall that for z/ = ^ ^ a = 1 and for A = 1 both models contain the ACDM model 
as a special case. The traditional Chaplygin gas with a = 1 is based in higher-dimensional 
theories 521]. The generalized Chaplygin gas with < a < 1 has been related to a Born- 
Infeld type approach 13]. The time dependence of the pressure is described by 



P 



G p 
G+^p 



P 



p 
P 



p 
P 



Gp 

Gp + ^ 



(6) 
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With 



P = -Q{p + p) 



e 



7 = 1 + ^ 
P 



(7) 



the adiabatic sound speed can be written as 

p p 
P P 



1 

- + 1^ 
2 



For the energy density we have 

3 



P 



a 



2 

l-2i/ 



H 



StiG 



1 



(9) 



The deceleration parameter q 



-1 — takes the form 

B /'ao\|(l-2i/) 



1 _ _B f'ao'i 

V a / 



Its present value go is 



qo = - 



1 - ^ 



s 1 + * 



1 + 1 



(10) 



1 - 2qo 

The value Oqcc at which the transition from decelerated to accelerated expansion occurs, is 
given by 

2 

^ \ 3(1-2^) 



g = <^ 



ao 



ao 



2AJ 



(12) 



Denoting the redshift parameter at which the acceleration sets in by Zacc, provides us with 



1 + Za 



ao 



1 - 2go^ 

1 + go 



- 1 . 



In terms of go the Hubble function in becomes 

1 - 2go + 2 (1 + go) 



H 



X \ 1-21- 

3 



while the corresponding energy density is {3Hq = SnGpo) 



P_ 

Po 



l'\ 1-21^ 

9 



l-2go + 2(l + go)(^^ 



2 



For the equation of state parameter 2 we obtain 



P 



l-2go 



1 - 2go + 2 (1 + go) (^) 



00^5(1-2'^) 



(13) 



(14) 



(15) 



(16) 



which implies 



7 = 1 + 



P_ 2(l + go)(^) 



P l-2go + 2(l + go)(^^)^^'-''^) 
The same relations hold for the generalized Chaplygin gas with the replacement 1 — 2v 

2(1 + a). 



(17) 
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III. PERTURBATIONS 



For a pressure p oc —p^Q the corresponding linear perturbations, denoted by the hat 
symbol, are 

18 



P 



e p 



p. 



Quantities without a hat are background quantities. The perturbations fllSp are non- 
adiabatic. Namely, 



p 




7^0. 



(19) 



Adiabatic perturbations are characterized hy p = ?p. It is the difference from p = jp which 
makes the perturbations non-adiabatic. The expression ([2]) for the pressure coincides with 
an equation of state p = p{p) oc — p'^+^Z^ only in the background. On the perturbative level, 
Eq. ([2]) cannot be reduced to p = p{p). Use of the relations allows us to write 



p 



p p 



p + p p p + p 



p\p e 



or, with the abbreviations 



P 



p 



p + p 



-D = 3H 
P P \ P 




(20) 

(21) 
(22) 



Both the combinations P — on the left-hand side and ^ — § on the right-hand side of 



(122|) are gauge-invariant, while the quantities P, D, p and 6 by themselves are not gauge- 
invariant. Obviously, the basic quantities for the study of the non-adiabatic perturbation 
dynamics are the energy-density perturbation p and the perturbation G of the expansion 
scalar. This suggests starting with the first-order energy conservation equation, while the 
perturbation of the expansion scalar is governed by the first-order Raychaudhuri equation. 
The general line element for scalar perturbations is 



ds^ = - (1 + 20) dt^ + 2a^F„dtdx° + [(1 - 2^/;) + 2E„^] dx"dx^ . 
The perturbed 4-velocity is described by 



and 



= Uq 



(23) 



(24) 



(25) 
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which defines the velocity perturbation v. A choice v = corresponds to the comoving 
gauge. It is also useful to introduce 

X = a'(E-F) . (26) 



The combination f + x is gauge-invariant. It is convenient to describe the perturbation 
dynamics in terms of gauge-invariant quantities which represent perturbations on comoving 
(superscript c) hyp ersurf aces. These are defined as 

- = - + v , -^ = ^+v , ^ = ^+v . (27) 
P p 6 6 p p 

In our case we have 



p ®j^yp p 0+i/£ 



(2^ 



Recall that a constant bulk-viscosity coefficient corresponds to z/ = 0. The perturbed energy 
balance may be written 

+m(^-l^\+l-[Av + Ax) = Q, (29) 
^P + P J \P + P PP + Pj 

where A is the three-dimensional Laplacian. From the momentum balance we have in first 
order 

^ + f + ?) + = . (30) 



p+p p+P 

In terms of the quantities introduced in (1271) . the balances (129|) and (!30|) may be combined 
into 

^ -3/7^-^ + 6^ = . (31) 



P + Pj P P + P 



With 



a more compact form of 0311) is 



D^^^, (32) 
P + P 



D"" -3H-D' + = . (33) 
P 



The expansion scalar 6 is governed by the Raychaudhuri equation 

6 + ^6^ + 2 - cj2) - u'^^-A + 4TrG{p + 3p) = 0. (34) 

Up to first order, the perturbed Raychaudhuri equation can be written in the form 

6^ + -66^ + -0^^' = , (35) 
3 a'' 6 
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where 



p 

P + P 



(36) 



It is through the Raychaudhuri equation that the pressure gradient comes into play. The 
formulation of the perturbation dynamics in terms of is particularly appropriate in the 
present case, since via (cf. fl28|) and (l36l) ) 



p 
P 



QC 



(37) 



the perturbation ©'^ of the expansion scalar is directly related to the pressure perturbation. 
Use of fl37|) in Eq. fl33|) provides us with a direct relation between the pressure perturbations 
and the energy-density perturbations, 



P 



,2p V P 



{31 



The pressure perturbation consists of a term which is proportional to the energy-density 
perturbations D'^, but additionally of a term proportional to the time derivative D'^ of V^. 
Pressure perturbations are not just proportional to the energy-density perturbations as in 
the adiabatic case. There is an additional dependence on the time derivative of the energy- 
density perturbations. The relation between pressure perturbations P'^ and energy pertur- 
bations D'^ is no longer simply algebraic, equivalent to a (given) sound-speed parameter as 
a factor relating the two. The relation between them becomes part of the dynamics. In a 
sense, is no longer a "local" function of D'^ but it is a function of the derivative D'^ as 
well [53|. This is equivalent to p = p{p,p). It is only for the background pressure that the 
familiar dependence p = p{p) is retained. 

Combining Eqs. (jHSD and (155]) and transforming to the /c-space, we obtain (using 

the same symbols as in the coordinate space) the second-order equation 

1 p k"^ 



+ 3H 



9H' 



2 p 

3 2p 9 7p H^a^ 



3\2 p^ ' 



P 



'^Yp^' + ^'^-lw^a-^p \2p 



= . 
(39) 



It is obvious, that the pressure perturbations give rise to contributions both in the brackets 
that multiply D'^ and D'^. For comparison, we also write down the corresponding equation 



9 



for the generalized Chaplygin gas (subscript c): 



- 9H' 



2 p 
- + a- 
.3 P. 



--a(l + a)7---- 1-3- 
D p b p \ p 



a 



p k'^ 



p m'^a? 



(40) 



IV. NUMERICAL IMPLEMENTATION AND RESULTS 



For the numerical implementation it is convenient to use 

6, = P- = ^D'^ , 
P 



(41) 



instead of D'^. The subscript v stands for viscous and was introduced to distinguish the 
perturbations in our viscous dark-fluid model from those of the corresponding generalized 
Chaplygin gas. In terms of the scale factor a, Eq. fl39l) then takes the form 



(42) 



where a prime denotes a derivative with respect to a and the coefficients /v(a) and gv{a) are 



/v (a) = - 
a 



- -6- + 3z/- - 

2 p p 3 7p _H^a^ 



and 



9v [a = — ^ 



3 15 p 9 p 



P 



2 p 2 p 



P 



lp2 



p\ k"^ 



IP 



p J H^a? 



(43) 

(44) 
M and (IE 



respectively. The quantities if, ^ and 7 as functions of a are given in ([T 
respectively. The present value of the scale factor is set to = 1 in the the numerical 
calculations. 

The perturbation equation for the generalized Chaplygin gas that corresponds to ( l42l) is 



(5;' + /c(a)5; + (7c(a)5c = 0, 



(45) 



with 



and 



/c(a) 



3 15p p 

3a- 

2 2 p 



{a) = — ^ 



p 2p2 



P 



(46) 



(47) 



respectively. With a = —{u + |) all the terms in ( l46l) and f j471) that are not multiplied 
by /c^ coincide with the corresponding terms in (l43i) and (jSj), respectively. The quantities 
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if, ^ and 7 as functions of a are also given by ( fT4l) . ( |T6l) and ( fTTI) . respectively, with the 
replacement 1 — 2z/ = 2(1 + a). The quantities (H3l) . (I44l) . (H6l) and ( 1471) are functions of the 
scale factor and depend on the parameters k, u, go, Hq. We shall restrict ourselves to z/ < |. 
We recall that sound propagation in the viscous fluid model is governed by a combination 
of the k"^ terms both in (143!) and in (jHj). In the Chaplygin-gas model, on the other hand, 
the sound speed square is given by — a^, the factor that multiplies the k'^ term in (j47j) . In 
contrast to the viscous fluid case there does not appear a k'^ term in (H^ . Eq. (H5i) with (H^ 



and (H7|) reproduce the basic perturbation equations in [15| and |21 |. 

At early times, i.e., for small scale factors a <C 1, both the equations (HTj) and (H5|) coincide 
and take the asymptotic form 

5" + A5'-A5 = 0, (a«l) (48) 

for all parameters go, ^ and for all scales. Here, 6 may either be 6^ or 6c- The solutions of 
are 

6{a < 1) = cia + C2a~^/^ , (49) 

where ci and C2 are integration constants. This means, at early times, the two models 
are indistinguishable. In particular, the non-adiabatic contributions to the viscous model 
are subdominant on all scales. Moreover, for a ^ 1 we can also consider our model to 
be indistinguishable from the ACDM model. This will allow us to follow the evolution of 
all models from the same initial conditions. We shall use the fact that the matter power 
spectrum for the ACDM model is well fltted by the BBKS transfer function |54!]. Integrating 
the ACDM model back from today to a distant past, say z = 1.000, we obtain the shape 
of the transfer function at that moment. The spectrum determined in this way is then 
used as initial condition for our viscous model. This procedure is similar to that described 
in more detail in references 



55 



56|. The numerical integration of Eq. (H2l) and Eq. (j35 



was performed with the help of the mentioned ACDM initial conditions where we used the 



parameters of the WMAP5 and 2dFGRS best-flt data sets 571]. 

In Figs. [1] and [2] the density fluctuations for the viscous model are compared with those 
of the generalized Chaplygin-gas model for different values of the relevant parameters. Al- 
though identical in the background, both models are qualitatively very different at the per- 
turbative level. The density perturbations in the bulk-viscous scenario are well behaved at 
all times, while there appear instabilities or oscillations in the GCG model, depending on the 
parameter a. The latter behavior was the main reason for discarding these models, except. 
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possibly, for very small values of a. Fig. [T] shows the blow up of the Chaplygin-gas density for 
i/ = 0(a = — 1/2). This reproduces a result of |26|. For z/ = — l(« = l/2), the Chaplygin- 

. Neither 



gas model predicts (unobserved) oscillations (cf. Fig. [2]), as was also found in 
of these unwanted properties hold for our viscous model. This coincides with the results of 
jsol . Both models coincide for early times, confirming our previous analytical result, that 
non-adiabatic contributions are negligible in the past, but become relevant at a later period. 
The non-adiabatic contributions are essential to avoid the mentioned unrealistic features of 
Chaplygin-gas models. 

The results for the matter power spectrum are shown and compared with the 2dFRGS 
and SDSS samples for different parameters v and go in Figures [3ll9l Two main features are 
observed here: (i) the bulk-viscous model is different from the ACDM model for all parameter 
choices. In particular, this is also true for the case ^ = ~\ (corresponding to a = 0) which 
has a background equation of state p oc —p. (ii) For a certain range of the parameters v and 
go, the model is in agreement with both the 2dFGRS and SDSS data samples. There occur 
neither oscillations nor instabilities. Negative values of v are generally preferred. The more 
negative v is, the more negative values of go are compatible with the data. Large negative 
values of v correspond to large positive values of a. This is consistent with the results of 



17 



18| . which 



recent studies on the perturbative behavior of generalized Chaplygin gases 
prefer large values of a as well, although, on the other hand, pure Chaplygin-gas models 
suffer from causality problems for a > 1. 

Finally we perform a x^-analysis and compare the values for the bulk-viscous model with 
the corresponding numbers for the ACDM model. The results are summarized in Tab. [B Our 
model turns out to be competitive with the ACDM model for go ^ —0.1, with a minimum 
of around go ~ 0. It is expedient to point out, however, that the result of the comparison 
depends strongly on the priors and does not seem to be a good indicator for the quality of 
the model. 



V. SUMMARY 

Unified models of the dark sector of the Universe may well be compatible with current 
observational data if this sector behaves as a bulk viscous fluid with a bulk-viscosity coeffi- 
cient C, oc . In the homogeneous and isotropic background they coincide with generalized 
Chaplygin-gas models with a = — ^ — z/, which are known to provide an adequate description 
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FIG. 1: Absolute values (logarithmic scale) of density fluctuations as function of the scale factor 
a for = (a = —1/2) and qq = —0.5 for different scales. The values of k are k = 0.5 (top left), 
k = 0.7 (top right), k = 1 (bottom left) and k = 1.5 (bottom right), all in units of hMpc^^. Solid 
curves represent the bulk viscous model, dashed curves the corresponding GCG model. Notice that 
both models are always different, except at very early times. 
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FIG. 2: Absolute values (logarithmic scale) of density fluctuations as function of the scale factor 
a for = — 1 (a = 1/2) and 50 = —0.5 for different scales. The values of k are k = 0.5 (top left), 
k = 0.7 (top right), k = 1 (bottom left) and k = 1.5 (bottom right), all in units of hMpc^^. Solid 
curves represent the bulk viscous model, dashed curves the corresponding GCG model. Notice that 
both models are always different, except at very early times. 
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FIG. 3: Density power spectrum for the bulk-viscous model with v = 0.25 (solid curves) and the 
ACDM model (dashed curves). From top to bottom the curves represent cases with go = —0.4, 
qo = —0.2, go = and go = 0.1. The curves are compared with 2dFGRS data (top) and and SDSS 
data (bottom). 





FIG. 4: Density power spectrum for the bulk-viscous model with = (solid curves) and the 
ACDM model (dashed curves). From top to bottom the curves represent cases with qq = —0.4, 
qo = —0.2, qo = and go = 0.1. The curves are compared with 2dFGRS data (top) and and SDSS 
data (bottom). 





FIG. 5: Density power spectrum for the bulk-viscous model with v = —0.25 (solid curves) and the 
ACDM model (dashed curves). From top to bottom the curves represent cases with go = —0.4, 
qo = —0.2, go = and go = 0.1. The curves are compared with 2dFGRS data (top) and and SDSS 
data (bottom). 




FIG. 6: Density power spectrum for the bulk-viscous model with v = —0.5 (solid curves) and the 
ACDM model (dashed curves). From top to bottom the curves represent cases with go = —0.4, 
g^o = —0.2, go = and go = 0.1. The curves are compared with 2dFGRS data (top) and and SDSS 
data (bottom). 




FIG. 7: Density power spectrum for the bulk-viscous model with v = —1.5 (solid curves) and the 
ACDM model (dashed curves). From top to bottom the curves represent cases with go = —0.4, 
g^o = —0.2, go = and go = 0.1. The curves are compared with 2dFGRS data (top) and and SDSS 
data (bottom). 




-2.0 -1.8 -1.6 -1.4 -1.2 -1.0 -0.8 



Logiok/h(Mpc-^) 



4.8 




-2.0 -1.8 -1.6 -1.4 -1.2 -1.0 -0.8 -0.6 



Logiok/h(Mpc-i) 

FIG. 8: Density power spectrum for the bulk-viscous model with v = —3 (solid curves) and the 
ACDM model (dashed curves). From top to bottom the curves represent cases with = —0.4, 
go = —0.2, go = and go = 0.1. The curves are compared with 2dFGRS data (top) and and SDSS 
data (bottom). 
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FIG. 9: Density power spectrum for the bulk-viscous model with v = —5 (solid curves) and the 
ACDM model (dashed curves). From top to bottom the curves represent cases with = —0.4, 
go = —0.2, go = and go = 0.1. The curves are compared with 2dFGRS data (top) and and SDSS 
data (bottom). 



TABLE I: Comparison of the x^'Values for the bulk-viscous model and the ACDM model for 
different parameters v and gg- 



V 


90 


(2dFGRS) 


(SDSS) 


0.25 


-0.3 


2830.33 


3776.76 




-0.2 


351.59 


459.76 




-0.1 


75.07 


72.89 







39.17 


54.79 




0.1 


35.40 


72.43 




0.5 


37.12 


98.76 





-0.3 


982.51 


2225.7 




-0.2 


238.79 


413.22 




-0.1 


85.97 


100.08 







47.90 


51.42 




0.1 


37.40 


56.36 




0.5 


37.12 


98.76 


-0.25 


-0.3 


570.22 


1453.14 




-0.2 


183.90 


331.18 




-0.1 


81.20 


98.02 







48.92 


52.39 




0.1 


38.44 


53.44 




0.5 


37.12 


98.76 


-0.5 


-0.3 


388.61 


997.17 




-0.2 


147.49 


256.43 




-0.1 


75.76 


87.27 







47.84 


51.74 




0.1 


38.57 


53.08 




0.5 


37.12 


98.76 


-1.5 


-0.3 


150.60 


299.40 




-0.2 


80.13 


104.61 




-0.1 


52.88 


56.75 







41.54 


50.06 




0.1 


36.96 


56.69 




0.5 


37.12 


98.76 


-3 


-0.3 


74.85 


95.75 




-0.2 


51.57 


55.53 




-0.1 


41.64 


49.91 







37.34 


55.27 




0.1 


35.69 


64.08 




0.5 


37.12 


98.76 



of the SNIa data [20|. While Chaplygin-gas type cosmologies show a pathological behavior 
on the perturbative level, the corresponding bulk viscous model is well behaved and its mat- 
ter power spectrum is compatible with the 2dFGRS and the SDSS observational data. It is 
the non-adiabatic character of the viscous fluid perturbation dynamics which is responsible 
for this difference. In a sense, our viscous model can be seen as a non-adiabatic general- 
ized Chaplygin gas. The model is observationally distinguishable from the ACDM model. 
Large negative values of v are preferred {v ^ —3). For certain parameter combinations a 
X^-analysis favors our model over the ACDM model. However, we found a discrimination on 
this basis not sufficiently convincing since it depends strongly on the priors. Our findings 



confirm and improve previous results on bulk- viscous cosmological models [49| . But we con 
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sider the present study as preliminary, since it does not explicitly take into account a baryon 
component. A corresponding generalization is currently under investigation. 
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